import control
import matplotlib.pyplot as plt

# 定义分子分母系数
Ka = 1999
num = 50 * Ka
den = [1, 60, 1100, 6000, 50 * Ka]

# 创建传递函数对象
sys = control.TransferFunction(num, den)

# 计算阶跃响应
t, y = control.step_response(sys)

# 绘制阶跃响应曲线
plt.plot(t, y)
plt.xlabel('Time (seconds)')
plt.ylabel('Response')
plt.title('Step Response of the System')
plt.grid(True)
plt.show()